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While band insulators are usually described in wavevector space in terms of fully filled bands, they are 
sometimes also described in terms of a complementary Wannier picture in which electrons occupy localized, 
atom-like orbitals. Under what conditions does the latter picture break down? The presence of irremovable 
quantum entanglement between different sites can obstruct a localized orbital description, which occurs in 
systems like Chern and topological insulators. We collectively refer to such states as Quantum Band Insulators 
(QBIs). Here we report the theoretical discovery of a filling-enforced QBI - that is, a free electron insulator 
in which the band filling is smaller than the minimum number dictated by the atomic picture. Consequently 
such insulators have no representation in terms of filling localized orbitals and must be QBIs. This is shown 
to occur in models of certain cubic crystals with non-symmorphic space groups. Like topological insulators, 
filling-enforced QBIs require spin-orbit coupling. However, in contrast, they do not typically exhibit protected 
surface states. Instead their nontrivial nature is revealed by studying the quantum entanglement of their ground 
state wavefunction. 


I. INTRODUCTION 

An early triumph of quantum mechanics was the realiza¬ 
tion that the metallic or insulating behavior of crystalline sys¬ 
tems could be understood in terms of the partial or full fill¬ 
ing of energy bands in wavevector space. For a filled band, a 
complementary, real space viewpoint has electrons localized 
at points in the unit cell, which in conjunction with the nuclear 
cores effectively lead to atomic picture of the insulating state. 
This was formalized by the work of Wannier et al. 1 which 
described a procedure to go from fully filled bands to a col¬ 
lection of symmetric and exponentially localized real-space 
functions, Wannier functions, akin to a collection of isolated 
atomic orbitals. Whenever such a description is applicable, 
the system can be understood as an atomic insulator (AI). In 
order to preserve crystalline symmetries, the centers of these 
orbitals occupy ‘Wyckoff positions’ in the unit cell and trans¬ 
form into each other under the space group (SG). 

However, not all band insulators are AI, meaning they do 
not admit such a real space depiction - we refer to such in¬ 
sulators as ‘quantum band insulators.’ The best known ex¬ 
ample is the Chern insulator in two spatial dimensions 2,3 for 
which there is an obstruction to constructing localized Wan¬ 
nier functions 4,5 . The obstruction stems from their nontriv¬ 
ial band topology - a collection of localized atomic orbitals 
cannot produce emergent chiral edge states. Similarly, topo¬ 
logical insulators (TIs) present an obstruction to realizing lo¬ 
calized Wannier states that preserve time reversal (TR) sym¬ 
metry 6 . Here we report the theoretical discovery of a ‘filling- 
enforced’ Quantum Band Insulator (feQBI) that shares several 
features in common with 3D TIs 7 - including (i) the absence 
of localized Wannier functions that preserve symmetry, and 
hence no real space picture of the insulating state and (ii) real¬ 
ization in a free electron Hamiltonian in which spin-orbit cou¬ 
pling is a crucial ingredient. In contrast to TIs, however, the 
feQBI is required by the electron filling v, i.e., the number 
of electrons per unit cell. A TI can generally be tuned, pre¬ 
serving symmetries but with a bulk gap closing, into a trivial 
AI at the same electron filling by modifying the Hamiltonian, 


while for the feQBI the gap, symmetries, and electron density 
alone forbid a competing AI state. Specifically, the feQBI is 
realized at a filling smaller than twice the minimum number 
of Wyckoff positions required to respect all SG symmetries. 
SG symmetries play a key role for feQBIs, akin to topological 
crystalline insulators (TCI) 8 . For example, we will show that 
in SG No. 199 with spin-orbit coupling (SOC), a feQBI oc¬ 
curs at a filling of 4 electrons per unit cell. However, based on 
the atomic picture the minimum AI filling is 8. In fact, feQBIs 
can only occur in 4 out of the 230 space groups (No. 199, 214, 
220 and 230) all of which are non-symmorphic, cubic crystals. 
Since cubic symmetry is always broken at a surface, feQBIs 
do not necessarily have protected surface states - their nontriv¬ 
ial character will be revealed instead via a suitably engineered 
entanglement-based approach. 

At first sight, the feQBI may seem related to a Mott insula¬ 
tor, which also occurs at ‘fractional’ filling. However, we em¬ 
phasize the feQBI is a state of noninteracting electrons, unlike 
Mott insulators. Furthermore, general interacting arguments 
which determine allowed fillings for band insulators, like that 
of Hastings-Oshikawa-Lieb-Schultz-Mattis 9-11 and their gen¬ 
eralizations to non-symmorphic crystals 12,13 with spin-orbit 
coupling 14 , correctly predict 4 as the minimum allowed fill¬ 
ing for a band insulator with space group No. 199. 


II. ATOMIC VERSUS BAND INSULATORS 

In this work we focus on systems of non-interacting spin- 
1/2 electrons with particle number conservation, TR and SG 
symmetries. SOC plays an important role, so there is no spin- 
rotational invariance. For brevity, we set all lattice constants 
to 1. 

Given an SG Q, symmetry-related points in space form a 
crystallographic orbit, and the collection of connected orbits 
is called a ‘Wyckoff position’ 15 . Each Wyckoff position 
is labeled by a ‘Wyckoff letter’ (w =a,b,...), ordered by de¬ 
creasing degree of site symmetry. A crystallographic orbit in 
has ‘|W£|’ points within each unit cell. Wyckoff posi- 
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tions with higher site-symmetry have smaller |W£|. 

SGs are classified as being either symmorphic or non- 
symmorphic (NS). An SG is symmorphic if one can pick an 
origin ‘o’ for which all symmetry elements can be decom¬ 
posed into a point group operation which leaves ‘o’ invariant 
followed by a lattice translation. By construction, the crys¬ 
tallographic orbit of ‘o’ is the Bravais lattice, so |W a | = 1. 
By definition such an origin is absent in a NS SG, and hence 
|Hf | > 2 if £ is NS. 

Wyckoff positions are intimately related to AIs, which are 
real space specification of a system. In the strict AI limit, 
electrons are localized to a set of points, though the points 
need not coincide with the real atomic sites in the problem. 
Nonetheless, to respect SG symmetries the set of points must 
correspond to a collection of crystallographic orbits, each in 
some Wyckoff position, i.e. the set of points constitute the 
lattice points of a SG-symmetric lattice. Since TR has to be 
respected locally, the system can be an AI only when the num¬ 
ber of electrons per point is even. Note that in a ‘valence bond’ 
state, the electron positions can be smoothly deformed to the 
center of mass of each valence bond, which reduces to the 
above point-like picture. AIs are therefore only possible at the 
‘atomic fillings’ v Al = 2 \^w I f° r some set of Wyckoff 

positions {w}. As Wf always has maximal site symmetry, it 
is always (one of) the position(s) with smallest |W£|. Hence 
the minimum filling for an AI is = 21 Wf |. For symmor¬ 
phic crystals, z / A1 n = 2, while for NS crystals z / A1 n > 4. 

Free electrons in a periodic potential are also convention¬ 
ally studied in momentum space via electronic band struc¬ 
tures. SGs constrain the properties of the bands since they are 
required to furnish little group irreducible representations (ir- 
reps) at high symmetry momenta. It is well known that bands 
can be forced to become degenerate due to the dimension¬ 
ality of the irreps, or due to time-reversal pairing at TR in¬ 
variant momenta 16 . In addition, the irreps at two end points 
of a high-symmetry line are not necessarily compatible, and 
such ‘compatibility relations’ can give further constraints on 
the possibility of isolating a set of bands. 

The minimal filling for forming a band insulator, ^f^ d , is 
identical to the smallest number of connected bands that can 
be isolated from all others. In particular, the class of free elec¬ 
tron problems contains those with no hopping between sites. 
The ground state of such a Hamiltonian is a strict AI, and 
hence z /^f n n d < z/ A1 n . In fact, for all SGs one can show that 
equality holds for spinless, TR-symmetric electrons 17 . 

Relaxing the condition on spin-rotation invariance, can the 
minimal filling for a band insulator be strictly smaller than that 
for an AI, i.e. can < z/ A1 n ? Here we answer this question 
in the affirmative: there exist SGs for which TR-symmetric, 
SOC band insulators can appear at a filling smaller than any 
AI. A band insulator realized at such sub-atomic filling is nec¬ 
essarily non-atomic, and hence the name feQBI. 


III. EXAMPLE OF A FEQBI 

The feQBI examples were found by a systematic study of 
all the 230 SGs. A particularly interesting observation is that 


for almost all SGs, the number of sites required to furnish 
any Wyckoff positions is always an integer multiple of | Wf |, 
which implies AIs are only possible at fillings z^ AI = n z/^ n for 
some integer n. There are, however, four exceptions: No. 199, 
214, 220 and 230 (Table I). We term such SGs as ‘Wyckoff- 
mismatched’. 

The simplest of the four Wyckoff-mismatched SG is 
No. 199, which has |W a 199 | = 4 and therefore ^ A1 n = 8. Yet, 
by studying the spinful little group irreps (Appendix A) and 
compatibility relations, we discovered a band insulator can be 
formed at a filling v = 4, i.e. ^f” d = 4 < z/^ n and hence 
proving the existence of feQBI. As a concrete example, we 
provide here a simple tight-binding model for this. We con¬ 
sider a set of sites in W b 199 , which contains six sites in each 
primitive unit cell. On the l-th (l = 1, 2,..., 6) site of the 
cell at x, we consider an s-orbital with electron operator 
where s =t, i corresponds to up and down spin quantized 
along the crystalline z-axis. The positions of the sites within 
each unit cell r l can be generated by the SG elements, start¬ 
ing from r 1 = (1/8,1,1/4). The tight-binding Hamiltonian 
is given by 

H =/^'t (t 5 S ' S + iA {a z ) s , s ) fi s + h.c. (1) 
+ (symmetry-related terms), 

where ‘h.c.’ denotes Hermitian conjugate, and repeated in¬ 
dices are summed over. In particular, H is TR-symmetric 
when the parameters A, t are real, and the ‘symmetry-related 
terms’ include all terms generated under SG symmetries. Note 
that in writing Eq. (1) we have made a particular choice of the 
unit cell, and a more explicit form of the Hamiltonian and 
the symmetry transformation of the sites are provided in Ap¬ 
pendix A. 

As shown in Fig. lc, when X/t = 1/4, the system is a 
band insulator at filling v = 4. This simple observation alone 
establishes existence of feQBIs. 

Before proceeding we note the following points. Without 
SOC (A = 0), the system at v = 4 is semi-metallic (Fig. Id). 
This is in fact the manifestation of a more general result: no 
feQBI is allowed when spin-rotation invariance is restored 17 . 
Interestingly, the lower four bands are also completely flat, 
since Eq. (1) describes a nearest-neighbor tight-binding model 
defined on the frustrated hyper-kagome lattice. In contrast 
to the flat bands of, say, the kagome lattice, here they can¬ 
not correspond to localized, SG-symmetric orbitals due to the 
non-atomic filling. 


IV. ENTANGLEMENT SPECTRUM OF THE 
‘SPIN-ORBITAL ENTANGLEMENT’ CUT 

We have established that an insulator at filling v = 4 in No. 
199 cannot be deformed to an AI without breaking any sym¬ 
metries (even when gap-closing is allowed), and by contrast¬ 
ing with the case of spin-rotation invariant systems one sees 
that SOC is essential for feQBI. How does SOC enable such 
band insulators at non-atomic fillings? Here, we argue that 
the feQBI has irremovable entanglement between two spin 
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FIG. 1. feQBI example with space group No. 199. (a) Hyper- 
kagome lattice with nearest-neighbour bonds considered in Eq. (1). 
Spatial symmetry-respecting spin quantization axes can be assigned 
on each of the sites, which we indicate by colored arrows, (b) Irre¬ 
ducible Brillouin zone of the space group, for which any momentum 
in the first Brillouin zone is symmetry-related to a point inside. Mo¬ 
menta are given in units of 7r, setting the lattice constants to 1. (c,d) 
Example band structures for the tight-binding model in Eq. (1). E F 
denotes the Fermi level at filling v — 4. With spin-orbit coupling 
in (c), the system is insulating and forms a feQBI. When spin-orbit¬ 
coupling is switched off in (d), the lowest four bands are completely 
flat, but they touch the upper bands at H and the system becomes 
semi-metallic. 


‘species’, grouped in a way that each spin species is closed 
under spatial symmetries but interchange under TR. Such en¬ 
tanglement is irremovable as long as the filling is fixed and 
symmetries are respected, and hence serves as a complemen¬ 
tary view on the symmetry obstruction to defining any AI at 
the same filling. 

A standard technique to probe entanglement between de¬ 
grees of freedom (dof) is via an entanglement cut and its cor¬ 
responding entanglement spectrum. An entanglement cut is a 
partition of the system into two subsystems A and A, corre¬ 
sponding to Hilbert spaces Ha and Ha such that the Hilbert 
space of the full system ‘factorizes’: H = Ha ® Ha- Most 
familiarly, A and A are spatial regions, though other partitions 
are possible 18-22 . If A is entangled with A, then the reduced 
density matrix pA = Tr^(|\F) will be mixed. The spec¬ 
trum of pA , often expressed via the ‘entanglement Hamilto¬ 
nian’ pa = e~ HE , contains universal fingerprints of quantum 
order. A TI has an intriguing gapless entanglement spectrum 
at a spatial cut. 23 In contrast, because an AI can be adiabat- 
ically deformed into an un-entangled state by strongly local¬ 
izing the orbitals, its ‘entanglement gap’ can be deformed to 
infinity while keeping all symmetries. 

The SG symmetries protecting feQBI, however, are broken 
by any naive spatial cuts. To make maximal use of the symme¬ 
tries, we introduce a ‘spin-orbital entanglement cut’ (SE cut) 
that both respect all symmetries and directly probe the irre¬ 
movable entanglement between spin-species. The full Hilbert 
space is a product over the Fock spaces H^ la of f x i a , where 
x labels spatial location, l labels orbitals and a = ff, Jj- labels 
the spin states, i.e. H = 0 xZcr ^xZa- Here we subsystem 
A be the collection of spin ‘up’ electrons instead of a spatial 


region, and the definition of A can be further engineered to re¬ 
spect all the symmetries at hand. The entanglement spectrum 
of the SE cut can be found by tracing over, say, the occupancy 
states |n) x ^. 

We would like all the symmetries ‘ g\ including the SG, 
to generate symmetries U g for the entanglement Hamiltonian 
H S e, i.e. U^HseU s = H$e (Appendix D). Since electron spin 
also transforms under spatial symmetries, one must choose the 
spin quantization axes in a symmetric fashion. Whether such 
choices exist or not depend on both the SG and the lattice re¬ 
alization, and we present a detailed discussion on this in Ap¬ 
pendix D. Here, we simply note that such SG-symmetric SE 
cut is possible for SG No. 199 (Fig. la). In the following we 
let A = ft denote the set of dof in the spin-up sector defined 
by the SE cut, and A = JJ, its complement (analogous to ‘left’ 
and ‘right’ in a typical spatial entanglement cut). 

TR exchanges ff, J|, so its role here is analogous to spatial 
inversion for spatial entanglement cut 23 . As detailed in Ap¬ 
pendix D, TR is realized as a unitary symmetry on H S e' while 
it is anti-unitary, the exchange of the two sub-systems intro¬ 
duces an additional complex conjugation in its proper defini¬ 
tion. The particle number N is also a good quantum number 
of H$e , and the sum of charges of the two TR-paired Schmidt 
states is fixed by the total electron number. Consequently TR 
is manifested as a particle-hole like symmetry. 

With the symmetries of H$e in hand, we now specialize 
our discussion to free electrons. The entanglement Hamilto¬ 
nian Hse of a free system is also free 18,23-25 , and can be re¬ 
duced to the single-particle entanglement Hamiltonian h $e = 
log(C _1 — 1), where Cij = (\F|is the correlation 
matrix and and i,j are restricted to dof in ff, and hence run 
over all of space. Note that C is nothing but the projector onto 
occupied bands further projected onto the ff dof. In contrast to 
a spatial cut, which introduces a boundary, Hse behaves like a 
local Hamiltonian throughout the entire bulk. 

By construction /zse respects all the SG symmetries. In par¬ 
ticular, it is periodic and hence one can discuss the entangle¬ 
ment energy bands ^(k), where l labels the bands. In ad¬ 
dition, (k) are continuous since the physical ground state 
corresponds to a band insulator. At physical filling z/, it is 
straightforward to verify that only the lowest v of the entan¬ 
glement bands can have (k) < +oo, so in what follows 
we consider only the lowest bands 1 < l < v. Hse has a 
gapped unique ground state if and only if h $e is gapped about 
the entanglement Fermi level ^| E = 0. 

For a strict AI, TR is realized locally so /zsE(k) consists of 
exactly v j 2 bands with £sE(k) at — oc, and all other bands are 
at +oo. In particular, the many-body entanglement ground 
state, constructed by filling all the z/se = v j 2 bands below 
£| e , is gapped and unique. The infinite entanglement gap is 
consistent with the intuition that there is no spin-orbit entan¬ 
glement in this limit. Away from the strict atomic limit, h se 
always could become gapless at d? E , even for an AI, but the 
question is whether Hse must be gapless. In the following we 
argue that the v = 4 feQBI symmetric under SG No. 199 nec¬ 
essarily has a gapless Hse- 

We argue by contradiction. Assume Hse is gapped with a 
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unique ground state. Because TR acts as a particle-hole sym¬ 
metry, the filling of the entanglement ground state must be 
half the physical filling, u S e = zV2 = 2, implying bands 1, 
2 lie isolated below the Fermi level: (k) 

— ^SE (k) < 

0 < (k) < 4E(k). However, we find that at k = P = 

( 7 r, 7 r, 7r) the four bands must transform under a ID and 3D 
irreps. It follows from the fact that the collection of irreps at 
P carried by the v entanglement bands must be the same as 
those of the v physical bands. This is because h$E is unitar- 
ily related to C , the occupied band projector further projected 
onto the f space; the latter projection respects the SG symme¬ 
tries and so leaves irreps unchanged. The feQBI here has one 
ID irrep and one 3D irrep, and this is a robust property of the 
system: the counting 4=1+3 cannot be altered without clos¬ 
ing the gap, and is independent of any particular tight-binding 
model used. This forces ^ 2 )(P) = which is incom¬ 

patible with an entanglement gap and hence a contradiction. 



FIG. 2. Spin-orbital entanglement cut and entanglement spectrum, 
(a) Schematic of the entanglement cut, in which the occupancy of one 
spin species (say Jf) is traced over, (b) Example spin-orbital entan¬ 
glement cut spectrum for the feQBI model of space group No. 199. 
The entanglement band structure inherits little group representations 
from the physical bands, and time-reversal symmetry is realized in a 
‘particle-hole’ manner. Note that at TR invariant momenta, different 
irreps of the physical bands can be forced to be degenerate due to TR 
pairing. Such degeneracies can be lifted in the entanglement bands, 
as seen at T and N of the example we gave. These together force an 
unavoidable gaplessness about e F (dashed line). 


V. IMPLICATIONS AND FUTURE DIRECTIONS 

While we have focused on SG No. 199 to illustrate the ex¬ 
istence and novelty of feQBIs, it is natural to ask if feQBI 
can exist in systems with other SG symmetries. In a recent 
work 14 we derived a general lower bound on the minimal fill¬ 
ing for all SGs, but the bounds reported there are not neces¬ 
sarily tight. For free electron problems, tight bounds can be 
found by symmetry analysis of possible band structures. We 
have carried out a comprehensive investigation on the possi¬ 
bility of hosting feQBIs for all SGs, and the results will be 
presented elsewhere 17 . Here we simply note that, in fact, with 
SOC a^in d = ^min f° r b ut ^e f our Wyckoff-mismatched 
SGs: No. 199, 214, 220 and 230. These four SGs admit fe¬ 
QBI at zz^ d = 4,4, 8 and 8 respectively. Interestingly, it ap¬ 
pears that feQBI arises from the interplay between Wyckoff- 
mismatch and the spin-1/2 nature of electrons. Understanding 
this link is an interesting problem we leave for future investi¬ 


gation. 


TABLE I. feQBIs in the four Wyckoff-mismatched space groups. 
‘Wyckoff multiplicities’ denotes the number of lattice points per 
primitive unit cell required to form a space-group symmetric lattice 
corresponding to one of the Wyckoff positions. a Al and zy Band respec¬ 
tively denote the electron fillings consistent with atomic and band 
insulators, and their discrepancy corresponds to feQBIs. 


Space group 

Wyckoff multiplicities 

u A1 

^Band 

No. 199 (72i3) 

4,6,12 

4N\{4} 

4N 

No. 214 (/4i32) 

4,4,6,6,8,12,12,12,24 

4N\{4} 

4N 

No. 220 (/43d) 

6,6,8,12,24 

4N\{4,8,20} 4N\{4} 

No. 230 (IaSd) 

8,8,12,12,16,24,24,48 

8N\{8} 

8N 


Conventional wisdom holds that one can study symmetry 
protected topological phases by finding obstructions to adi- 
abatically deforming the system to an un-entangled product 
state (like the AI) while respecting all symmetries. There were 
known examples of a ‘filling-enforced’ obstruction to such a 
deformation if one insists on using a particular tight-binding 
realization of the SG 26 . These examples, however, are thought 
to be trivialized once symmetry-allowed modifications to the 
lattice models are included, equivalent to bringing down dof 
from infinity. Our result indicates there are filling-enforced 
obstructions to such deformation (adiabatic or not) given only 
the symmetries and filling, independent of the lattice or con¬ 
tinuum realization. While the topological class of the feQBI 
may well fall under established classifications (say K-theory, 
which allows for band addition, and field theories, which do 
not incorporate the filling 7,27 ), the ‘filling-enforced’ nature of 
the feQBI appears to be new, and it would be worthwhile to 
revisit these programs with ‘filling’ in mind. 

We close by discussing possible material realization of fe¬ 
QBI. As indicated in this work, one should consider a system 
with significant SOC in SG No. 199 or 214 with v = 4 (dis¬ 
regarding core electrons tightly bounded to atoms), or in SG 
No. 220 or 230 with v = 8. In addition, the energy bands rele¬ 
vant near the Fermi energy should contain the correct symme¬ 
try irreps needed for realizing feQBI (Appendix A). Although 
we do not yet have a concrete realization in mind, here we 
point out the hyper-kagome iridate NasIraOg comes close to 
providing the ingredients needed: in a hypothetical structure 
promoting the simple cubic bravais lattice to body-centered 
cubic (by altering atom positions), the relevant TB Hamil¬ 
tonian can be adiabatically connected to the model given in 
Eq. (1) in certain parameter region (Appendix C). 


VI. CONCLUSION 

In conclusion, a SOC electronic system with TR and SG 
symmetries can form a band insulator at a filling that is in¬ 
compatible with any AI. This also defies the usual perspective 
that a TI or TCI can be understood as a ‘twisted’ version of 
an AI, having been rendered non-atomic due to certain pat- 
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terns of band inversion at high-symmetry momenta 7 . More 
generally, our result implies that for certain symmetries there 
are electron fillings for which a symmetry protected topo¬ 
logical phase is the only option if the system is symmetric, 
gapped and short-ranged entangled. We also point out ma¬ 
terials exhibiting the hyper-kagome structure are potentially 
experimental playground to realize feQBIs. 

The spin and orbital dof are intrinsically entangled in a fe- 
QBI. This was exposed through a ‘spin-orbital entanglement 
cut’, which can be engineered to make maximal use of all 
symmetries and is also applicable to interacting many-body 
systems. Curiously, a given feQBI system can simultaneously 
be a conventional TI (for instance, the model in Eq. (1) is 
a strong TI), but this is not a necessary consequence of the 
filling-enforced nature of feQBIs (Appendix B). It is worth¬ 
while to point out, however, that the Wyckoff-mismatched 


SGs all have SG symmetries that are broken by the existence 
of any surface, so finding physical bulk signatures of feQBIs 
is an important open question. 
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Appendix A: Symmetries of space group No. 199 

In this appendix we provide supplementary information on 
the symmetry properties of SG No. 199, an explicit tabulation 
of the feQBI tight-binding example discussed in Eq. (1), and 
also an additional tight-binding example. 

1. Symmetries of SG No. 199 

Up to lattice translations, No. 199 contains 12 symmetry 
elements. We parameterize the SG element g by g = {R g = 
p e -i0 g -L, | ^ where / = +1 (—1) for proper (improper) 
rotation. While these can be easily determined from the in¬ 
formation listed in ITC 15 , we tabulate them in Table II for 
completeness (‘el’, denotes the symmetry operation indexed 
by the same number in ITC). 


SG No. 199 has three Wyckoff positions labeled by a,b and 
c, which can be realized by having 4, 6 and 12 sites respec¬ 
tively in each primitive unit cell. To aid visualization of W a 199 
and W b 199 , we present in Fig 3 representative lattice realiza¬ 
tions of them viewed from different angles. 

2. SG symmetric spin-texture 

Sites in W a 199 are invariant under a three-fold rotation, but 
the rotation axis is site-dependent. A SG symmetric spin- 
texture can be formed by polarizing spins along the corre¬ 
sponding axis on each site. Sites in W b 199 are invariant up to 
lattice translation under one of the two-fold screws, and so the 
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TABLE II. List of symmetry elements in Q/T for SG No. 199, pa¬ 
rameterized by g = {R g — /e _z6>9 ' L | t g}. 


el. 

I 

0g 

*9 

(i) 

+ 

(0,0,0) 

(0,0,0) 

(2) 

+ 

tt(0, 0,1) 

(i0,2) 

(3) 

+ 

7r(0, 1,0) 

(o, \) 

(4) 

+ 

7r(l,0,0) 

(!> !>o) 

(5) 

+ 

f 7s(M,i) 

(0,0,0) 

(6) 

+ 


(i §>0) 

(7) 

+ 

-1.-1) 

a, 0 , 2 ) 

(8) 

+ 

¥ts(-W,i) 

(0, \) 

(9) 

+ 

f 0M,1) 

(0,0,0) 

(10) 

+ 

f 7i(W,-l) 

(0,U) 

(11) 

+ 

— —(—1 -1 11 

(!> !>o) 

(12) 

+ 

— —(—1 1 —11 

3 ’ ’ ’ 

a, 0 , 2 ) 



FIG. 3. (a,b) View of representative sites forming W a 199 along high 
symmetry axes. There are 4 sites in each primitive unit cell and hence 
8 sites in the conventional cell. (c,d) Same as (a,b) but for W b 199 . 


spin-texture has the spin polarized along the corresponding 
screw axis. The spin-quantization axes are uniquely defined 
for sites in W a 199 and W b 199 (up to a global choice of up vs 
down), but to be self-contained we list them in Table III. 

We also clarify here a subtle point (though non-essential for 
our arguments) : if the screws were ‘intrinsically’ NS, no site 
should be taken back to itself nor its lattice-translation images. 
This, however, is not the case for SG No. 199 as the screws are 
not intrinsic: given any of the screws, one can choose an ori¬ 
gin such that the ‘screw’ factorizes into a point-group rotation 
followed by lattice translation. No. 199 is nonetheless non- 
symmorphic because no common origin can be picked to ren¬ 
der all space group elements symmorphic. SGs like No. 199, 
which are NS but do not contain any intrinsically NS element, 


TABLE III. Spin-quantization axes corresponding to the SG sym¬ 
metric spin-texture. The sites are labeled by l in the same or¬ 
der as they are listed in ITC 15 . The axes are parameterized as 
fii — (sin 6i cos (pi, sin 6i sin <j>i, cos 6i), and we denote #[in] = 
cos -1 (l/\/3). 


l 

Wi" 

0i 

4>l 

l 

wr 

0i <f>i 

1 

°[ 111] 

7t/4 

1 

7T/2 

0 

2 

0[111} 

—37t/4 

2 

7r/2 

7T 

3 

7T - 0[111] 

37t/4 

3 

7T/2 

7t/2 

4 

7T — 0[in] 

—7t/4 

4 

7t/2 • 

—7r/2 




5 

0 

0 




6 

7r 

0 


are known as ‘exceptional’ NS SGs, and in fact there are only 
two of them out of the 230 SGs: No. 24 and No. 199. 


3. Little group irreducible representations 

In Table IV we list the symmetry eigenvalues of the dif¬ 
ferent irreducible representations of the little group at high 
symmetry momenta T, H and P. The set of isolated four 
bands forming feQBI corresponds to p^ 0 at T and 

Pp lc ^ 0 Pp 3a ^ at P. Irreps at other high symmetry momenta are 
fixed once these are specified. 

Appendix B: feQBI tight-binding examples 
1. Model given in Eq. (1) 

One can deduce the transformation of tight-binding (TB) 
dof by assigning them to some actual locations in real space. 
In Eq. (1), we assume the sites have coordinates in W b 199 at 
x = 1/8, with a single 8-orbital on each site. Note that the 
TB model is chosen such that it contains the irreps needed to 
realize feQBI, but feQBI is possible in any TB models that 
contains the irreps required. 

We label the sites (1 - 6) in the same order as they appear 
in ITC 15 (left - right). In writing their coordinates as r = 
x + with x a lattice vector, r l depends on the choice of 
the unit cell. Here, we take the primitive lattice vectors as 
ai = (1/2,1/2,1/2), a 2 = (0,1,0) and a 3 = (0,0,1), and 
choose r l such that r l are contained inside the parallelepiped 
defined by the three primitive lattice vectors. With this choice 
of unit cell, the coordinates are given by 

r 1 = (1/8,1,1/4); r 2 = (3/8,1,3/4); 

r 3 = (1/4,9/8,1); r 4 = (1/4,7/8,1/2); (Bl) 

r 5 = (0,1/4,1/8); r 6 = (0,3/4,3/8). 

The transformation of the TB basis is parameterized by writ¬ 
ing g(r l ) = ni&i + r l , where {rii} is a triplet of integers. 
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TABLE IV. Symmetry eigenvalues of the irreducible little group rep- TABLE V. Transformation of tight-binding sites, labeled by l, under 

resentations at high-symmetry momenta. We denote uj = e -12 ^ 3 . the symmetry elements (el.), n denotes — n. 


(a) T = (0,0, 0) and H = (2tt, 0,0) : 


el. 

P^ 

P? b) 

P? c) 

(i) 

{1,1} 

{1,1} 

{1,1} 

(2) 

{*, -*} 

{*,-*} 

{*,-*} 

(3) 

{*, -*} 

{*,-*} 

{*,-*} 

(4) 

{i, -i} 

{*,-*} 

{*,-*} 

(5) 

{-(V, -UJ*} 

{—1, —a;} {-1,-cc*} 

(6) 

{-(V, -UJ*} 

{—1, —uj} { —1, —u>*} 

(7) 

{-cc,-cH} 

{-!.-«} {-!,-"*} 

(8) 

{-UJ, -UJ*} 

{-1,-w} {-1,-w*} 

(9) 


{!,«*} 

{l,cc} 

(10) 

{w,w*} 

{!."*} 

{l.w} 

(11) 

{w,w*} 

{!."*} 

{l.w} 

(12) 


{!,«*} 

{l,cc} 


el.\ l 
( 1 ) 
( 2 ) 

(3) 

(4) 

(5) 

(6) 

(7) 

( 8 ) 

(9) 

( 10 ) 
( 11 ) 
( 12 ) 


1 2 3 4 5 6 

(1; 000) (2; 000) (3; 000) (4; 000) (5; 000) (6; 000) 

(2; 020) (1; 021) (4; 021) (3; 020) (5; llO) (6; 120) 

(2; 110) (1; 110) (3; Ill) (4; 110) (6; 000) (5; 010) 

(1; 121) (2; 122) (4; 122) (3; 122) (6; III) (5; III) 

(3; OlO) (4; llO) (5;2l0) (6; llO) (1; OlO) (2;0l0) 
(4; 112) (3; 223) (6; 323) (5; 212) (1; HI) (2; H2) 

(4; Oil) (3; III) (5; 102) (6; Oil) (2; OlO) (1; Oil) 

(3; lOl) (4; 210) (6; 210) (5; IlO) (2; 100) (1; 100) 

(5; 2ll) (6;2II) (1; 211) (2; HI) (3; Oil) (4; HI) 

(6; 211) (5; 221) (2; 321) (1; 211) (3; 100) (4; 211) 

(6; 322) (5; 322) (1;332) (2; 222) (4; HI) (3; 222) 

(5; 101) (6; III) (2;2Il) (1; III) (4; OlO) (3; IlO) 


TABLE VI. A full list of terms in the feQBI tight-binding example 
given in the main text. The terms are parameterized as in Eq. (B2). 


el. 

(la) 

PP 

p^ 

(1 c) 

pp 


p^ 3a) 

k • Sx 

l' 

A 

A 

l 

(i) 

1 

l 

1 


{1,1,1} 

0 

4 

t 

(0,0, A) 

1 

(2) 

-1 

-l 

-1 


{-1,1,1} 

0 

3 

t 

(0,0, A) 

2 

(3) 

-1 

-l 

-1 


{-1,1,1} 

0 

4 

t 

(0,0,-A) 

2 

(4) 

-1 

-l 

-1 


{-1,1,1} 

-kz 

3 

t 

(0,0,-A) 

1 

(5) 

-1 

—UJ 

—UJ* 

{- 

-1, -w, -w*} 

2 {P>x T Py T kz') 

6 

t 

(A, 0,0) 

3 

(6) 

-1 

—UJ 

—UJ* 

{- 

-OJ*} 

\(k x + ky + kz) 

5 

t 

(A, 0,0) 

4 

(7) 

-1 

—UJ 

* 

—UJ 

{- 

-1, — UJ, — UJ*} 

0 

6 

t 

(—A, 0,0) 

4 

(8) 

-1 

—UJ 

* 

—UJ 

{- 

-1, -w, -w*} 

ky + k z 

5 

t 

(—A, 0,0) 

3 

(9) 

1 

* 

UJ 

UJ 


{1,UJ,UJ*} 

— \ (kx + ky + kz) 

2 

t 

(0, A, 0) 

5 

(10) 

-1 

—UJ* 

—UJ 

{- 

-1, —uj, — UJ*} 

0 

1 

t 

(0, A, 0) 

6 

(11) 

-1 

* 

—UJ 

—UJ 

{- 

-1, — UJ, — CJ*} 

~ \ {kx + ky + k z ) 

2 

t 

o' 

o' 

6 

(12) 

-1 

* 

—UJ 

—UJ 

{- 

-1, -w*} 

— ky 

1 

t 

o 

o' 

5 


While this is readily computable, in Table V we tabulate the 
transformation of site l as (/'; n \, n- 2 , v/ 3 ) under the symmetry 
elements. 


To construct a Hamiltonian H that is fully symmetric un¬ 
der SG Q, one can start from a single (non-symmetric) term 
H 0 , like the one given in Eq. (1), and take the summation 
H = J2 g eg dHog- 1 ■ For reader’s convenience, we give an 
explicit form of H in the momentum space here. To fix nota¬ 
tion, we specify a generic term in a periodic Hamiltonian for 
SOC electrons as follows: 

=e -*-SxfW(t Ssla + i(A • a) s , s )fl s + h.c. 

(B2) 

where Sx is also a lattice vector. In Table VI, we specify the 
Hamiltonian by providing a list of all terms in this notation. 


2. An alternative ‘minimal’ model 

Since = 4 for SG No. 199, the ‘minimal’ tight- 

binding model (corresponding to local Hamiltonians) should 
have at least eight bands. This is indeed possible by consider¬ 
ing a pair of orbitals on each of the sites furnishing W a 199 . 
However, to match the representation content required for 
forming a feQBI, the TB dof have to transform in a more in¬ 
trigue manner under the SG symmetries. 

To this end, note again that each site in W a 199 is invariant 
under a three-fold rotation along a site-dependent axis, and 
hence it can pick-up an orbital phase of 1 , uj = e ~ l27T ^ 3 or c o 2 
under the three-fold rotations. A spin-1/2 polarized along the 
three-fold rotation axis will also pick up a phase of e _z27r / 6 
for spin-l| and e l27T / 6 for spin-J|. In particular, we consider a 
TB model formed by the {uj, ff) and {uj 2 , J|) orbitals on each 
of the symmetry-related sites belonging to W a 199 . Note that 
the two on-site orbitals also interchange under TR. 

We again parametrize the TB model in a similar manner as 
in Eq. (B2), but here we use an unconventional basis for the 
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Pauli matrices (corresponding to indices s and s'), namely we 
write ( t jm ) = ( fjlt t fm ) U(6i,<j)i), where 

and are respectively the creation operators for the (cc, -ft) 
and (cc 2 , ft) orbitals on the Z-th site in the unit cell with co¬ 
ordinate x, and U(Qi,(j)i ) is the site-dependent unitary trans¬ 
formation relating the crystalline z-axis to the SG symmetric 
spin-texture (Table III), i.e. had we picked the ft and ft states 
from the same orbital then would just correspond to the 
creation operator written in the basis of crystalline z-axis. 

We consider a TB model with four sites in each unit cell 
starting with r 1 = (0, 0, 0). Similar to the previous example 
the coordinates for l = 1 ,..., 4 can then be determined 
from the SG symmetries (same choice of unit cell as before), 
with the sites labeled in the same order as in ITC. The terms in 
the eight-band TB Hamiltonian are tabulated in Table VII. In 
particular, TR invariance of the Hamiltonian still requires A, A 
to be real, although the A term is no longer spin-independent 
in this choice of basis. One can check that the system forms 
a feQBI when, say, t' 2 / t'i = 2. The corresponding electronic 
and entanglement bands are plotted in Fig. 4. 


TABLE VII. An alternative eight-band feQBI tight-binding example. 


k • £x 


r 

t 

A 

l 

0 


4 

t'i 

(0,0,0) 

1 

0 


3 

A 

(0,0,0) 

2 

0 


2 

0 

(0,40) 

1 

ky 


4 

0 

(0,40) 

3 

2 (J^X ky 

k z ) 

3 

0 

(0,0, -4 

1 

|(-fcc +ky ~ 

- k z ) 

4 

0 

(0,0, -4 

2 

-k z 


4 

0 

(0,0,4 

1 

-k z 


3 

0 

(0,0,4 

2 

~ky 


2 

0 

(—40,0) 

1 

0 


4 

0 

( 4 o,o) 

3 

2 ( k X ky 

-k z ) 

3 

t'2 

(0,0,0) 

1 

\(k x —f ky 

k z ) 

4 

t'2 

(0,0,0) 

2 


3. Hyper-kagome lattice and feQBI in the 
Wyckoff-mismatched space groups No. 214,220 and 230 

In this subsection we show that the simple TB model in 
Eq. (1) can be used to establish the existence of feQBIs in 
all the four Wyckoff-mismatched SGs: No. 199, 214, 220 and 
230. We reproduce in Table I the multiplicities of the Wyckoff 
positions for these four SGs. Note that not all | | are integer 

multiples of | Wf | for these SGs. 

To this end, we first focus on SG No. 214. As listed in Table 
I, |W 214 | = 4 and so = 8 for SG No. 214. Now we note 
that the coordinates listed in Eq. (Bl) are chosen to coincide 
with W 214 . Keeping only nearest neighbor bonds, the system 
is in a hyper-kagome structure (three-dimensional network of 
corner sharing triangles) and each site is four-coordinated. 



FIG. 4. Alternative feQBI example for SG No. 199. Note that the 
symmetry representation content is identical to the example given in 
the main text, (a) Band structure with t ' 2 / 1 [ = 2. Note that the 
model has an extra ‘sub-lattice’ symmetry (i.e. there exists a unitary 
U such that UH\JJ — —HE) and hence the bands are symmetric 
about — 0, but this is purely an artifact of the simple TB model 
and is not essential, (b) Corresponding entanglement band structure. 


Demanding TR and spatial symmetries of SG No. 214, the 
symmetry allowed terms take the form 

H =/x S,f (t S s > s + i{ A x a z + A 2 ( Cf x - o v )) s ’ s ) f^ s (B3) 
+ h.c. + (symmetry-related terms). 

Eq. (1) of the main text corresponds to the special case when 
A 2 = 0, and so it can just as well be regarded as an example 
of an insulating state inSGNo. 214 with filling v = 4 < 8. 
This establishes the existence of feQBI for SG No. 214. For 
completeness, we also plot in Fig. 5 the phase of the system 
described by Eq. (B3) at filling v — 4. 

Next we consider the centrosymmetric SG No. 230, which 
can be viewed as SG No. 214 supplemented with spatial in¬ 
version. In particular, W 230 is the union of W 214 and W^ 14 
, which are inversion-copy of each other. Consequentially, 
one can simply take the TB model in Eq. (B3) together with 
its inversion copy to construct a TB model symmetric under 
SG No. 230 (defined on two inter-penetrating hyper-kagome 
lattices that are decoupled). Hence one can construct a TR 
symmetric insulator with SG No. 230 at v = 8. Since 
8 < 2|W 230 | = 16, this is also a feQBI. 

Finally, we note that SG No. 220 is a subgroup of No. 230. 
In particular, 8 < 2|>V 220 | = 12 and so the same model con¬ 
structed for No. 230 is also a feQBI example for No. 220. 


4. Strong and weak indices of feQBI examples 

The strong and weak indices were found by computing the 
Z 2 indices for the six TR-symmetric planes spanned by two of 
the reciprocal lattice vectors in the Brillouin zone containing 
either T or G^/2 (where G* is the remaining reciprocal lattice 
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FIG. 5. Plot of band gap A E for the hyper-kagome model in 
Eq. (B3) at filling v — 4. For each pixel of the plot, A E is deter¬ 
mined by sampling 30 x 10 3 momentum points along the high sym¬ 
metry lines and in the first Billouin zone respectively. ★ corresponds 
to the system for Fig. l(e,f) of the main text and A corresponds to 
that of Fig. 6. 



FIG. 6. Plot of surface band structure against the surface crystal 
momentum ky for the model in Eq. (B3). We take Xi/t = 0.5 and 
M/t — 0.15. E v denotes the bulk Fermi energy. Slabs with thick¬ 
ness of 20 surface-adapted unit cells are used and results for different 
surfaces, corresponding to different surface normals n, are shown. 
Note that with open boundary conditions, there are surface states but 
they do not traverse the bulk gap. (a) Periodic boundary condition 
for ft = z. (b) Open boundary condition for h — z. (c) Open 
boundary condition for ft ~ y + z. (d) Open boundary condition for 
fi ~ x + y + z. 


vector not spanning the plane). The strong index was found 
to be uq = 1 and 0 respectively for the SG No. 199 feQBI 
models in Eq. (1) and the one specified in Table VII, while 
both of them have weak indices characterized by the vector 
Gzv = \i'iGi = 2ttx. As such, these models actually feature 
surface states (on the appropriate surfaces), even though they 
are not required by the filling-enforced nature emphasized in 
this work. 

We note that, however, it is not necessary for a feQBI to si¬ 
multaneously process nontrivial Z 2 indices. An explicit exam¬ 
ple for a feQBI with trivial Z 2 indices is the inter-penetrating 
hyper-kagome model constructed for SG No. 220 and 230 
(cf. previous subsection). There the various Z 2 indices are 
simply given by twice of the corresponding indices of the 
model in Eq. (1), and therefore must be trivial. 

As another example, observe that Fig. 5 features two ‘is¬ 
lands’ of insulating phases, separated by a gapless region. 
This is not a coincidence. The island marked by ★ corresponds 
to the strong TI model given in the main text, and the one 
marked by A is in fact neither a strong nor weak TI. As such 
the mentioned gapless phase between the two islands should 
actually feature a nodal semi-metal. We plot in Fig. 6 the sur¬ 
face band structures of A on various surfaces, all of which are 
gapped. 


Appendix C: Hypothetical structure for spin-orbit-coupled 
hyper-kagome material Na 3 Ir 3 Os 

In this appendix we point out the observation that the exper¬ 
imentally synthesized hyper-kagome material Na3lr 3 0g (Na- 
338) is in close proximity to a feQBI phase, and can point to 
promising avenues for the experimental realization of feQBIs. 

Na-338 can be regarded as a hole-doped version of Na- 
438, a Mott-insulating hyper-kagome material well-known as 
a spin-liquid candidate 28 . In reality Na-338 crystallizes in SG 


No. 213, which is a primitive version of SG No. 214. The 
structure of the atoms are listed in Table VIII (adapted from 
Ref. 29). 

TABLE VIII. Measured structure 29 of Na 3 Ir 3 0 8 in SG #213 (Nota¬ 
tion as in ITC 15 ) 


Atom Wyckoff Free Parameter Representative 


Ir 

12 d 

y — —0.113 

(0.613,0.863,5/8) 

01 

8 c 

x = 0.114 

(0.114, 0.114, 0.114) 

02 

24e 

(x,y,z) 

(0.136, 0.907, 0.919) 

Nal 

4b 

- 

(7/8,7/8, 7/8) 

Na2 

8 c 

x = 0.257 

(0.257, 0.257, 0.257) 


Due to the strong SOC of Iridium and crystal field splitting 
coming from the local environment, the relevant states near 
the Fermi level can be (roughly) modeled after an effective 
J e ff = 1/2 Kramers pair centered at the Ir sites. This gives 
rise to 24 relevant energy bands around the Fermi level with 
an electron filling of v = 8. (Note that in this picture, we 
are regarding the other 48 bands arising from the J e ff = 3/2 
states of Ir as separated in energy and fully filled, contributing 
to ‘core electron states’.) 

While Na-338 is semi-metallic in reality 29 , here we con¬ 
sider a hypothetical structure of it in which the atom positions 
are modified to ‘enrich’ the spatial symmetries from SG No. 
213 (simple cubic) to No. 214 (body-centered cubic). In par¬ 
ticular, we change the free parameter associated with Ir in 
Table VIII from y = — 0.113 to y = 0, putting the Iridium 
atoms in W/ 14 . A possible assignment of the other atoms are 
provided in Table IX. 

In promoting the bravais lattice from simple cubic in SG 
No. 213 to body-centered cubic in No. 214, the electron fill- 
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TABLE IX. ‘Symmetry-enriched’ hypothetical structure of 
Na 3 Ir 3 0 8 in SG #214 


Atom Wyckoff Free Parameter Representative 


Ir 

12 c 

- 

(0,1/4,1/8) 

01 

8 a 

- 

(1/8,1/8,1/8) 

02 

24h 

y= 1/4 

(0,1/8,1/4) 

Na 

12 d 

- 

(0, 1/4, 5/8) 


ing per primitive unit cell is halved. Retaining only nearest 
neighbor bonds, the system can now be effectively described 
by Eq. (B3) with v = 4 and serves as a feQBI example if the 
parameters Ai and A 2 lie in the insulating phase indicated in 
Fig. 5. 

Note that in this discussion we only intend to give experi¬ 
mental context to feQBI, instead of proposing realistic mate¬ 
rial candidates. In particular, the ‘oxygen cage’ is significantly 
distorted from the ideal octahedron form in the structure tab¬ 
ulated in Table IX, and hence whether the effective spin-1/2 
picture still holds or not deserves scrutinization. Nonethe¬ 
less, we also note that the site-symmetry group for W 214 is 
the crystallographic point-group D 2 , which has only one 2- 
dimensional spinful representation. As such, as long as one 
can identify a Kramers doublet, well isolated in energy from 
other states, living on the sites of W 214 , the system is de¬ 
scribed by Eq. (B3) when restricted to nearest neighbor bonds. 


Appendix D: Discussions on the spin-orbital entanglement cut 

1. General discussion on time-reversal symmetry 

Consider an electronic system with particle number conser¬ 
vation. Let the total number of particles TV be even and |\P) be 
a TR symmetric many-body state. We fix the phase ambiguity 
such that T|\P) = |\P). Consider a basis {|i)} for sub-system 
ft and {|z)} for ft . In this basis, the ground state can be ex¬ 
panded as 

I*>=Z^jK>IJ> (Dl) 

ij 

Generally, M is a rectangular matrix. Since TR exchanges 
ft and ft, however, M is in particular square here. Singular- 
value-decomposition simply gives M = WXVY where W 
and V are unitary, and A is diagonal and positive semi-definite. 
The Schmidt states in this basis are given by 

l*> = E (IW°) (| 3)Vf a ) = E (D2) 

a a 

and in the same basis the reduced density matrices are p ^ = 
W 7 A 2 1-0 = MM* and p 4 = V*A 2 0 T = M T M*. 

Applying TR operator T, 

l*> = E A « (/»*) = E A « (I a f)b (l«r)t), 

a a 

(D3) 


and hence the two Schmidt states | a)^ and | af)^ are degen¬ 
erate in the entanglement spectrum. Note that while T itself is 
anti-unitary, the map relating the ‘ft’ Schmidt states with the 
‘ft’ Schmidt state is also anti-unitary. Altogether, TR becomes 
a unitary symmetry on the entanglement Hamiltonian. 

To see this more explicitly, let TR transform the basis by 

m = mb f\i) = mb < D4 ) 

where the unitary matrices U r and U T satisfy 

U f U f * = (-1)^, (D5) 

as required by T 2 = (— 1)^. TR invariance of T|^) = |^) 
implies 

M =U f M f (u f Y => p t = U f p\U f] . (D6) 

One can bring p ^ back to p ^ using a ‘spectral flattened’ ver¬ 
sion of M: M = WVY See that 

mVh-M = W*X 2 W T = p\, (D7) 

so altogether 

Pt = {U f M ] )p^{MU f] ), (D8) 

i.e. TR is realized as a unitary symmetry on the entanglement 
Hamiltonian, as claimed. 

In particular, since [TV, 7~] = 0, we also have 

N (|« t ) t ) = tN\a)t = (N — N a1t ) (|a t )^), (D9) 

and hence the charges of the paired states are symmetric about 
TV/2, i.e. TR is manifested in a ‘particle-hole’ manner. Note 
that since TR is now realized unitarily, a state with charge TV/2 
can be paired with itself under TR. 

To illustrate the ideas discussed above, we provide here a 
simple illustrative example. Consider a two site system with a 
particular TR symmetric state 

1*) = (ViVu+ mfUflx. - 7V1V2V) 

(DIO) 

+ s ( e flplf + e * flbb) 1 °)’ 

and under TR, -b and -b — Hence, 

= (a* fyi+- bfipb) i°) 

7 (Dll) 

+<5*(e/i t t / 2 t t + |0). 

Choosing the global phase such that T|\I/) = |^), we take 

rt, /T, G M. 

See that 

=Seftpt t |0>t; 

mini*) = (-a/iV - io) t ; 

v 7 (D12) 

1(01/24.1^} = ( — 1^7/it) l°)t; 

MfnfnW =**|o> t , 
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and so the reduced density matrix after tracing out occupancy 
states with label j. will have two degenerate eigenvalues of 
S 2 |e| 2 . These two Schmidt states, having charge 0 and 2 
respectively, form a TR pair. The remaining 2-dimensional 
block is (in the basis of /^|0)^ and |0)^): 


(P t ) 2 x2 = 


a 2 + P 2 \l\ 2 a/3 7 \ 

ah* P 2 \l\ 2 ) 


The corresponding Schmidt weights are 


(D13) 


Hamiltonian (even when SG-symmetric spin-texture cannot 
be defined). This is due to the existence of a ‘primary rota¬ 
tion axis’ on the sites: for these 27 point groups, a spin po¬ 
larized along the primary rotation axis is either flipped or left 
invariant (up to a phase) by any symmetry operations. A SG 
symmetry is then realized as an anti-unitary or unitary symme¬ 
try on the entanglement Hamiltonain depending on whether it 
flips the spin or not. In particular, we note that this is true for 
all the Wyckoff positions for the four Wyckoff-mismatched 
space groups. 


4 = Y+/3 2 |7l 2 ±a^hpw, < D14 > 

which are generally non-degenerate, and each state is paired 
with itself under TR. This is allowed since they have charge 
1 = 2/2. 

2. General discussion on space group symmetries 

The SE cut is defined to partition dof by their spin-label 
ff, J|, which are picked with respect to some quantization axes. 
In the presence of SOC, the spin states on different sites are 
also related by SG symmetries. Therefore, to ensure the SE 
cut is SG symmetric we must carefully choose the spin quan¬ 
tization axes in a site-dependent fashion. 

A sufficient condition for the SE cut to respect all SG sym¬ 
metries is the existence of SG-symmetric spin-texture, i.e. the 
spins can be polarized in a site-dependent fashion while pre¬ 
serving all the SG symmetries. Equivalently, we demand the 
existence of site-dependent spin-quantization axes for which 
any spatial symmetries either rotates the spin about the axis, 
or takes the site to some other site. 

Generally speaking, such SG-symmetric spin-textures may 
not exist for a given lattice realization of a space group: 
whether or not they can be defined depends on the site- 
symmetry of the sites in the lattice. Since a generic point 
in space, belonging to the general Wyckoff position, is never 
taken back to itself under any SG element, the only possible 
obstruction comes from high-symmetry points (correspond¬ 
ing to high-symmetry Wyckoff positions) where any choice 
of spin-polarization breaks some site symmetries. 

Nonetheless, as long as one is concerned about the space- 
group symmetries, but not the specific lattice realization of the 
system, any obstruction to defining a SG symmetric spin tex¬ 
ture is largely technical. Intuitively, in the continuum it should 
be possible to adiabatically ‘punch-out’ the high-symmetry 
points without changing the phase. More concretely, one can 
always ‘split up’ a high-symmetry orbital into a small set of 
orbitals in the vicinity of the high-symmetry point in order to 
avoid any obstruction. In other words, one can approximate 
any desired orbitals on a lattice site by effective ‘molecular 
orbitals’ on sites belonging to the lower-symmetry Wyckoff 
positions. 

We also note that if the site-symmetry group is one of the 
27 non-cubic crystallographic point groups, then all SG sym¬ 
metries can be regarded as symmetries of the entanglement 


3. Single-particle entanglement Hamiltonian 

In the following subsections we specialize our discussion 
to free electron problems. Due to Wick’s theorem, the ground 
state entanglement spectrum of a free fermion system parti¬ 
tioned into A and A is completely captured by the single¬ 
particle correlation matrix Cij = where |^) is 

the many-body ground state and i,j, G A 18,23-25 . Since the 
spin-orbital entanglement cut preserves spatial symmetries, 
and in particular translation invariance, we can focus on cor¬ 
relation functions evaluated at a fixed momentum k. This is 
related to the filled-band projector V k by Ck ;a /3 = (^k 
where a,/3 are now collective index for basis, orbitals and spin 
dof, and run from 1,..., 2 m . In a proper choice of basis, the 
first m entries of V k corresponds to the spin ‘up’ sector in the 
SE cut, and the m x m restricted correlation matrix is simply 



is similarly defined. The definition is, of course, basis 
independent - a basis transformation of 7\ is always accom¬ 
panied by a corresponding transformation of the rectangular 
matrix projector. Since ^sE(k) is completely determined by 
, we focus on the properties of C £ in the following. 

As C £ is a product of projectors, we have eig(C£) G [0,1]. 
A (many-body) Schmidt state |A) with charge N\ is formed 
by populating N\ single particle eigenstates |c$) of the full 
correlation Hamiltonian C ^ = © k ^k' The eigenvalues of 
the reduced density matrix of such a state is a product over the 
eigenvalues of the occupied (occ.) and unoccupied (unocc.) 
states 18,23-25 : 

A 2 = II Ci II ~ c i) (D16) 

i€ occ. j£ unocc. 

and hence the highest weight Schmidt states are formed by 
filling all states of C ^ with q > 0.5, i.e. in this language 
one should think of filling the bands of C ^ from above, not 
below, and the effective ‘chemical potential’ is 0.5. Note that 
the single-particle entanglement energies are related to q here 
by Si = log(cG 3 — 1). In particular, if there are states with 
c l = 0.5, they contribute to A in the same way whether they 
are filled or not, and corresponds to degeneracy in the entan¬ 
glement ground state. Note also due to particle-hole charac¬ 
ter of the TR symmetry on the entanglement spectrum, the 
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highest weight Schmidt state, if unique, must have a charge 
corresponding to half of the real filling. 


4. Symmetry properties of the single-particle entanglement 
Hamiltonian 

Recall that a single-particle Hamiltonian symmetric under 
a spatial symmetry g satisfies 

U{?H k U° =H g - I (k) (D17) 

where is a unitary. The notion of a spatially symmetric 
spin-texture implies that on each site, one can specify a spin 
quantization axis such that the two spin eigenstates are decou¬ 
pled under all spatial symmetries. This implies there exists a 
k-independent unitary matrix Ws such that 

WsU^Wl = ^ ^ j . (D18) 

For simplicity, in the following we assume we work with such 
a basis in mind such that we already have U£ = ® . 

Since the band projector satisfies 

= u£p g -i W U£, (D19) 

the restricted correlation matrix transforms as 

C» = (t/ k »») T C s % k| ([/»*)’ , (D 20 ) 

i.e. it transforms as a spin-polarized system with the same SG 
symmetries. 

It remains to show how the original TR symmetry is man¬ 
ifested in Note that for the original Hamiltonian, TR 
symmetry implies H k = and therefore V-k = 

basis of Eq.(D15), one simply has Ur = 

—ia y ( 8 ) 1 mxmi and therefore C £ = . Using the prop¬ 

erties of the projectors and the duality discussed in Refs. 18 
and 20, one sees that if G (0,1) is an eigenvalue of 
then 1 — is an eigenvalue of C^ k , i.e. TR is now manifested 
as a ‘particle-hole’ symmetry between ^sE(k) and ^se(— k) 
(Fig. 2b and 4b). More explicitly, suppose C^.ip = ^Ck, and 
let 

$ = ( 0 1 ) v£ ^ ^ • (D21) 

One can then verify = (1 — Ck)'?/'. Hence, as long as 
^ 7 ^ 0 , it is an eigenvector of with eigenvalue 1 — Ck- 
In addition, one can check that ^ = 0 implies Ck = 0 or 1, so 
any Ck G (0, 1) is paired with a c_k = 1 — Ck, i.e. translating 
this relation into the spectrum of ft-sE(k), TR is reflected as a 
particle-hole like symmetry with £ S E(k) = —£se(— k). 


